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ABSTRACT 

We utilize observations of sub-millimeter rotational transitions of CO from a Herschel Cycle 2 
open time program (“COPS”, PI: J. Green) to identify previously predicted turbulent dissipation by 
magnetohydrodynamic (MHD) shocks in molecular clouds. We find evidence of the shocks expected 
for dissipation of MHD turbulence in material not associated with any protostar. Two models fit 
about equally well: model 1 has a density of 10 3 cm -3 , a shock velocity of 3 km s _1 , and a magnetic 
field strength of 4 pG; model 2 has a density of 10 3 ' 5 cm' 3 , a shock velocity of 2 km s^ 1 , and a 
magnetic field strength of 8 /rG. Timescales for decay of turbulence in this region are comparable 
to crossing times. Transitions of CO up to J of 8, observed close to active sites of star formation, 
but not within outflows, can trace turbulent dissipation of shocks stirred by formation processes. 
Although the transitions are difficult to detect at individual positions, our Herschel -SPIRE survey of 
protostars provides a grid of spatially-distributed spectra within molecular clouds. We averaged all 
spatial positions away from known outflows near seven protostars. We find significant agreement with 
predictions of models of turbulent dissipation in slightly denser (10 3 ' 5 cm -3 ) material with a stronger 
magnetic field (24 fi G) than in the general molecular cloud. 


1. INTRODUCTION 


The dust in molecular clouds shields the interior from 
ultraviolet radiation, lowering the kinetic temperature to 
about 10 K. For this temperature, the isothermal sound 
speed is only about 0.2 km s _1 . In contrast, the observed 
linewidths of molecular transitions indicate much larger 
veloc ity dispersions. When the se facts first became evi¬ 
dent, [GoIdrel^&lKwaS ( 1974 ) suggested that the large 
linewidths were caused by overall collapse of the molec¬ 
ular cloud, rejecting the alternative of turbulent motions 
because such motions shoul d decay on a tim escale similar 
to a crossing time. iZuckerman fe Palmerl (19741) noted 
that allowing all molecular clouds to collapse at free- 
fall would produce a star for matio n rate far higher than 
that of the Milky Way. IZuckerman fc Evans: (1974 1) sug¬ 
gested that the linewidths were indeed caused by tur¬ 
bulence, and that magnetic fields could slow the dis¬ 
sipation if the motions were sub-Alfvenic. The idea 
that the linewidths are indicative of turbulent motions 


is wi despread dMvers fc Gammiel 119991: ISolomon et al 
1~98^: iLarsonl 1 1981 1) . but simulations (Stone et all 1199 


Smith et all 120001 ) indicate that the turbulence should 
decay via shocks. Even when magnetic fields are in¬ 
cluded, simulations indicate that turbulence decays, 
but through magnetohydrodynamic (MHD) or C-sh ocks 
(|Mullanlfl971I iDraind fl98ft iDraine fc Me <ed 1 1 9931) . or 
through coherent, small-scale vortices (Godard et all 
I2014D . Strong evidence for decay has been found 
in diffuse and transl uc ent clouds using chemical in¬ 
dicators (Godard et al.1 I2014D and H 2 emission lines 
(|Falgarone et all 12005( 1. The latter method has also in¬ 
dicated evidence for decay on the outskirts of the Tau¬ 
rus molecular cloud (Goldsmith et all 12010 ). The idea 
that turbulence decays in denser molecular clouds is also 
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widely accepted, but direct evidence for the shocks in 
opaque molecular clouds has proven elusive. The dis¬ 
covery of such evidence would resolve a long-standing 
question. 

The sources of this turbulence are still somewhat un¬ 
clear. One source for renewing the turbulence is out¬ 
flows from forming stars. The outflows are hypersonic, 
but they may couple energy to the surrounding gas via 
magnetic disturbances. Thus, we may find evidence for 
enhanced dissipation in the surroundings of outflows, but 
we must avoid the outflows to avoid confusion. 

Much of the energy deposited in the gas by shocks is 
radiated away, and a likely route is through CO and H 2 
rotational transitions, as well as through enhanced emis¬ 
sion by atoms and ions, such as OI and CII. Attempts to 
find signatures of turbulent dissipation in lines of H 2 , CII, 
and OI are frustrated by the fact that these lines could 
also be excited by low intensity photodissociation regions 
( PDR s) which cover the surface of all clouds (Li et al.1 
120021) . A renewed opportunity to detect the signature 
of the shocks has arisen due to the combination of new 
models, focusing on the mid-J CO lines and the ability 
of the Herschel instruments (both HIFI and SPIRE) to 
detect these lines toward regions that are not directly 
affected by outfl ows. 

Models from iPon et ahl (20121 ) indicate that turbulent 
shocks enhance mid- J CO emission over th at from qui¬ 
escent PDR models (Kaufman et al.l 11999 ). Observa¬ 
tional evidence for turbulent dissipation has been found 
in the Perseus Bl-East 5 region by IPon et aD ( 2014 ). 
They chose this region of the Perseus molecular cloud 
as one with no significant heating sources other than 
shocks, cosmic rays, and interstellar radiation fields 
(ISRF); there are no signs of embedded young stel¬ 
lar objects. Measurements of the CO J = 5 —> 4 and 
J — 6 —> 5 transitions show valu e s enhanced over PDR 
models by iKaufman et al.l (19991) . IPon et all (20141) con¬ 
cluded that the filling factor of shocks toward this po- 
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sition was low (0.15%), and that the timescale for tur¬ 
bulent dissipation was about a factor of 3 shorter than 
the flow crossing time dPon et al.ll2015bh . This is an im¬ 
portant discovery, but it applies only to a single region 
and required a combination of data from different in¬ 
struments. A confirmation with other techniques in dif¬ 
ferent regions would add confidence that turbulent dissi¬ 
pation has finally been detected. Similar results have 
recently been obtained for apparently starless clumps 
within IRDCs bv lPon et al.l ( 2015al ). 

In this work we present a sample of SPIRE spec¬ 
troscopy of fields in molecular clouds including one near 
no protostar and seven near protostars. Using portions 
of these fields unaffected by outflows or nearby proto- 
stellar radiation fields, we examine them for evidence of 
turbulent decay, as indicated by increased luminosity in 
mid-J CO lines. 

1.1. Object Selection 

The objects that were analyzed were taken from the 
COPS-SPIRE sample (Green et al. in prep.) of low lu¬ 
minosity embedded protostars. These observations were 
not optimized for this project as they tend to have power¬ 
ful outflows which produce emission in very high- J lines. 
To avoid outflows, we selected objects with weaker out¬ 
flows and with outflows well confined to certain regions, 
focusing on the regions away from the outflows. In ad¬ 
dition, we obtained one region within a filament of the 
Taurus cloud that did not include any embedded source, 
thus providing a sample of the molecular cloud without 
any star formation. We refer to this observation as the 
isolated molecular cloud. 

Table □ describes the seven regions with embedded 
sources that were selected for this study. L1455-IRS3, 
IRAS03301-3111, and TMR1 were chosen because 
they can be easily spatially resolved in the SPIRE 
field. TMC1 and TMC1A have compact outflows 
that allow for analysis of the surrounding molecu¬ 
lar cloud within the field of view. L1551-IRS5 has 
a large, but well collimated, bipolar outflow in the 
North-East/South-West direction. L1157 has distinct 
collimated bipolar outflows in the North-South direction. 


2. OBSERVATIONS AND REDUCTION 

Data sets were collected by the Herschel Sya ce Ob¬ 
servatory , a 3.5-m space telescope dPilbratt et al.ll2010i l. 
using the SPIRE (Spectral and Photometric Imaging 
Rece iver) Fourier Transform Spectrometer dGriffin et al.1 
12010( 1. The SPIRE short-wavelength (SSW) spectrum 
covers the 194-313 /jm range over an array of 37 pixels 
(two of which were not functioning) that project beams 
at 33" separation. The SPIRE long-wavelength (SLW) 
spectrum covers the 303-671 jjl m range over an array of 
19 pixels, and these beams are separated by 51". For 
this work, we restrict the analysis to CO lines within 
only the SLW data sets, which cover CO transitions up 
to J = 9 —> 8. We also examined the SSW pixels for 
higher-J lines, but most of our sources did not show 
evidence for CO in the SSW range, and the higher-J 
CO lines we could identify could not be resolved from 
the outflows of the forming stars. 


TABLE 1 
Source List 


Source 

Cloud 

Dist 

RA (J2000) 

Dec 

L1455-IRS3 

Per 

250 

03h28m00.4s 

+30d08m01.3s 

IRAS 03301+3111 

Per 

250 

03h33ml2.8s 

+31d21m24.2s 

L1551-IRS5 

Tau 

140 

04h31m34.1s 

+18d08m04.9s 

TMR 1 

Tau 

140 

04h39ml3.9s 

+25d53m20.6s 

TMC 1A 

Tau 

140 

04h39m35.0s 

+25d41m45.5s 

TMC 1 

Tau 

140 

04h41ml2.7s 

+25d46m35.9s 

L1157 

Core 

325 

20h39m06.3s 

+68d02ml6.0s 

Isolated MC 

Tau 

140 

04h39m53.0s 

+25d45m00.0s 


Note. — List of protostellar sources and the isolated molecular cloud 
position used in this sample, sorted by RA. Distances are measured in 
parsecs. 


We separate each data set by pixel, including the 
object itself and the area around it. The continuum 
in each pixel was first calculated using an automated 
routine written in Interactive Data Language (IDL) 
and was then subtracted from the original spectrum 
to create a flat (continuum-subtracted) spectrum. The 
flat spectra from each pixel in the combination were 
averaged together to improve the signal-to-noise ratio. 
The CO lines in each average were fitted by a Gaussian 
function which yields the intensity integrated over the 
line. This intensity was converted to the units used in 
the models (erg s -1 cm~ 2 arcsec -2 ). Upper limits for 
the non-detections were calculated by taking the noise 
level and multiplying by 3 to get the 3cr value. The CO 
J = 7 —> 6 line flux estimate is particularly uncertain as 
it is blended with the [C I] line which is very prominent 
in the isolated molecular cloud field. Figure |T| shows 
an example of the blended CO J = 7 —>■ 6 and [C I] 
3 p 2 —t 3 Pi lines and an overlay of the double-Gaussian 
fit that was applied to distinguish the fluxes from each 
line individually. 



Fig. 1.— Example of line fit for the blended CO J = 7 —> 6 and [C 
I] 3 P 2 —> 3 Pi lines from the combined average of the non-outflow 
pixels from Figure [6] 


2.1. Pixel Combinations 

For the isolated molecular cloud, we could average the 
spectra over all the pixels. For the embedded sources, 
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we used the following combinations of pixels: 

• Center - The central pixel (labeled SLWC3), con¬ 
taining the protostar. The SPIRE beam size is 
large enough that the continuum source is point¬ 
like and restricted to the center pixel. 

• Outflows - These pixels contain outflows driven by 
the protostar. For each observation, we compared 
the SPIRE pixels to ground-based CO J = 2 —> 1 
and J = 1 —> 0 contour maps from Kang et al. 
(submitted). An example of this overlay and com¬ 
parison can be found in Figure [2] 

• Non-outflows - This category includes all pixels not 
“contaminated” by the protostar or outflow emis¬ 
sion. These pixels should be the best representa¬ 
tion of the gas around a forming protostar, but not 
directly affected by the outflow. 



3. ANALYSIS 

We compare our data to t hree s ho ck m odels supplied 
by A. Pon (pers. comm.). iPon et akl (1201211 modeled in¬ 
tegrated intensities of various CO rotational transitions 
in units of erg s -1 cm~ 2 arcsec~ 2 , considering three val¬ 
ues for density: 10 2 5 cm -3 , 10 3 ° cm -3 , and 10 3 5 cm -3 . 
These values are based on the mo lecul ar cloud chem i- 
cal com position originally used in IKaufman fc Neufeldl 
(Il996al lblf. They considered shock velocities of 2 and 3 
km s _1 for a cloud temperature of 10 K, corresp onding 
to Mach 12 and Mach 17 respectively (IPon et al.1120121 ). 
These speeds are consistent with turbulence in molecular 
clouds but are significantly lower than the speeds of pro- 
tostellar outflows. Magnetic field strength (shown here 
in nG) in a molecular cloud is defined by 

B = bn k H 

where nn is the number density of hydrogen nuclei and b 
and k are fitting parameters. The value of k = 0.5 corre¬ 
sponds to a constant rat io of magnetic energy de nsity to 
thermal energy density (iMcKee fc Qstrikerll2007D and is 
the value assumed in these models. The strength of the 
magnetic field parallel to the shock front ranges from 3 
/xG to 24 /xG, characterized by b = 0.1 and b = 0.3. The 
shock thickness is smaller and the energy released in line 
radiation is larger for a weaker field. Therefore, these 
values were chosen as they are slightly l ower than, but 
still consistent with, those determined bv ICrutcher et all 

(l2M0|) . 

We also comp are our data to CO predictions from 
IPon et a kl (|2012f), which are based on the PDR models 
of IKaufman et akl (Il999f) . These PDR models have in¬ 
terstellar radiation field (ISRF) values specified in units 
of Habings (1.6 x 10 -3 erg cm -2 s _1 ). 

These models also assume line widths of 1.5 km s _1 , 
which are on the scale of the shock model velocities, and 
a total Ay of 10. They do not take into account CO 
freeze-out and also assume plane parallel geometry. The 
densities used for the PDR models correspond to the den¬ 
sities of the shock models and all assume a Go value 
of 0.31, which is reasonab ly consis tent with actual val¬ 
ues for Taurus. iFlagev et al.l (|2009D found a value of 0.2 
Habings. 



U (o'cuc) 

Fig. 2.— Example of outflow map from object L1157. Top: An 
overlay of the pixels of the SLW array from the Herschel telescope 
over the 12 CO J = 1 —> 0 contour map of the outflows from the 
object made by Kang et al. (submitted). The velocity ranges used 
for the blue and red components are —3.8 km s —1 to 4.0 km 1 
and 9.3 km s —1 to 10.6 km s~ 1 . respectively. Bottom: The SLW 
pixels overlaid on the 12 CO J = 2 —¥ 1 contour map. The velocity 
ranges used for the blue and red components are —0.6 km s —1 to 
4.3 km s —1 and 8.6 km s —1 to 9.4 km s —1 , respectively. 

We present these models in Figure [3] (isolated molec¬ 
ular cloud) and Figure [4] (grand average). The mod¬ 
els show lines ranging from J = 1 —^ 0 on the right to 
higher J levels toward the left. The model parameters 
are indicated in each figure panel. The plotted line fluxes 
are those taken from the averages obtained in the Cen¬ 
ter, Outflows, and Non-outflows composites: green (non¬ 
outflow pixels), purple (outflow pixels), and light blue 
(central protostar). For the isolated molecular cloud, 
all pixels are considered “non-outflow” (green). Upper 
limits to fluxes are indicated by an upside-down trian¬ 
gle of the appropriate color. The red points/curve show 
the intensities predicted for PDR models, while the dark 
blue points/curve are for the MHD shock model. These 
are clearly distinguished in the models in that the shock 
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spectra dominate over the contribution from PDRs for 
transitions from J = 5 —> 4 and higher (IPon et al.ll2012li . 




log Wovelength (/zm) 

Fig. 3. — Isolated Molecular Cloud: Observed CO emission lines 
from CO J = 8 —> 7 (left) to CO J = 4 —» 3 (right) from all pixels 
(green) plotted over two of the models they best fit. Triangles are 
estimated 3 a upper limits for non-detections. PDR models use a 
Go value of 0.31. Error bars are present but might be too small to 
distinguish. Top: Intermediate density (n = 10 30 cm -3 ), highest 
velocity (n = 3 km s —1 ) and lower magnetic field strength (B = 4 
pG). Bottom: Highest density (n = 10 3 ' 5 cm -3 ), lower velocity 
(v = 2 km s 1 ) and lower magnetic field strength (B = 8 //G). 

3.1. Isolated Molecular Cloud 

We begin with an isolated portion of the Taurus molec¬ 
ular cloud with no embedded protostar. It is centered on 
a filament, far from any apparent source. The emission 
is weak (we do not detect transitions above J = 6 —> 5) 
but the higher-J lines are clearly stronger relative to the 
lower-J lines than expected for a PDR model. Their 
strengths are reasonably consistent with both the 10 3 0 
cm -3 density, highest velocity (v = 3 km s _1 ), and lower 
magnetic field strength (B = 4 /zG) model and the high¬ 
est density (n = 10 3 ' 5 cm -3 ), lower velocity (v = 2 km 
s - 1 ) an d lower magnetic field strength (B = 8 /zG) model 
of IPon et aD (|2012D . as shown in Figure [3j Improved 
sensitivity to detect higher-J lines is required to further 
constrain the higher-J lines and constrain the choice of 
model. The fact that even the low Go PDR models over¬ 
estimate the observations of the J = 4 —> 3 line indicates 
that the PDR contribution to the higher J lines would 
be even less than shown in Figure 3. Even so, the ob- 



Fig. 4.— Grand average of all seven objects: Observed CO emis¬ 
sion lines from CO J = 8 —> 7 (left) to CO J = 4 —»• 3 (right) from 
the center (light blue), outflow (purple), and non-outflow (green) 
pixels plotted over the model they best fit to: highest density 
(n = 10 3 ' 5 cm -3 ), highest velocity (y = 3 km s —1 ), and higher 
magnetic field strength (B = 24 juG). PDR model uses a Go value 
of 0.31. Error bars are present but might be too small to distin¬ 
guish. 

servations exceed by an order of mag nitud e the PDR 
predictions for J = 6 —> 5. IPon et ah . (120141 ) considered 
other PDR models from different groups and found that 
they did not substantially change their conclusions that 
shocks were required to match the observations of the 
higher J lines. However, PDR models that produce sub¬ 
stantially higher excitation of the higher J CO lines could 
conceivably reproduce the observations. 

3.1.1. Dissipation Timescale 

The fact that the observations match the model for J 
values above J = 5 —> 4 allows us to estimate the dissi¬ 
pation timescale for the turbulence at the position of the 
isolated molecular cloud. We treat the column through 
the cloud defined by a single pixel of the SPIRE array as 
typical of any region in the cloud which does not contain 
a protostar. The turbulent energy in this column is given 

by 

E t = 0.5E m cr 2 D b d 2 (1) 

where E m is the mass surface density, a is the 3-D rms 
turbulent velocity dispersion, S2b is the beam solid angle 
in sr, and d is the distance to the cloud. The rate of 
energy dissipation through shocks is 

L sh = 47rd 2 F co fco (2) 

where Fco is the total observed line flux in all the CO 

lines and fco is the fraction of the total line emission 
emitted in CO. The flux in CO is the sum of all the 
intensities in the best fitting model times the solid angle, 
Hi,, (Fco = Ico^b)- The timescale becomes 

id = E t /L sh = fc ° E y a2 (3) 

where Ico is the sum of all the intensities of CO lines 

from the best fitting shock model in units o f erg s~ 3 
cm -2 sr -1 . The shock velocity in the models of lPon et all 
(2012 (1 is expressed in terms of the one-dimensional ve¬ 
locity dispersion (ctid) as v s h = 3.2<tid = (3.2/\/3)cr, so 
a = 0.54u s h For the two best fitting models, u s h = 2 or 
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3 km s _1 , n = 10 35 or 10 3 4 ° cm -3 , and Iqo = 2.11 or 
2.16 x 10“' erg s _1 cm -2 sr _1 . 

We determined £ m from maps o f exti nction base d on 
stellar reddening from iSchlaflv fo F inkbein ei) (120111) and 
far-infrared emission from ISchlegel et al.l (Il998h which 
yielded estimates of E(B — V) of 2.49 ± 0.02 mag 
and 2.89 ± 0.02 mag respectively, where the uncer¬ 
tainties clearly do not include systematics. We con¬ 
vert to Ay using Ry = 5.5, found to best mat ch gen¬ 
eral molecular cloud extinction laws dChapman et al.l 
12009H . Then we convert to £ m usin g the extinction effi¬ 
ciency from iWeingartner fc Drainel (120011) 0 for Ry = 
5.5, Case A (newer m odels). The results are 150 
pc~ 2 (ISchlafl v fe Finkbeinei! 12(1111) or 175 M 0 pc“ 2 
(ISchlegel et al.lll998lh If Ry = 5.5, Case B (older model) 
grains are used, the values are 206 and 240 M 0 pc -2 . For 
simplicity, we assume £ m = 150 M 0 pc“ 2 in what fol¬ 
lows, providing a lower limit to the dissipation timescale. 

For the model with u s h = 2, ta = 1.6 Myr; in the model 
with v s h = 3, <d = 3.3 Myr. Following the procedure in 
Pon et al. 2012 we u se the linewidth-size relation from 
[Solomon et al.l (|1987l ). r(pc) = (eriD/0.72km s -1 ) 2 where 
fj id = <t/ "\/3 to estimate a cloud radius of 0.77 pc for the 
v s h = 2 km s _1 model. Still following Pon et al. 2012 we 
estimate a flow-crossing time from t CT = 2R/ ctid of 2.3 
Myr for the u s h = 2 km s _1 model and 5.1 Myr for the 
v s h = 3 km s _1 model. The ratio of td to t CI is thus 0.94 
or 0.65 for the two models. Larger values of £ m would 
increase these ratios by a factor up to 1.6. This analysis 
shows that the turbulence dissipation timescale derived 
from observations and best-fit mod el is com parable to 
the flow-crossing timescale. IPon et akl (|2014H found a 
ratio of 1/3 (originally published as 3, but corrected in 
an Erratum) (IPon et al.ll2015bl) in the Perseus cloud. 

3.2. Embedded Sources 

The fields around the embedded sources can provide 
a sample of denser cores, where turbulence may be en¬ 
hanced. The challenge is to avoid mixing in positions 
toward outflows, which produce high velocity shocks and 
strong emission from high- J levels of CO. An example of 
a field around an embedded protostar (L1157) is shown 
in Figure [2] LI 157 provides a particularly strong test 
case, as the outflows (derived from CO J = 1 —» 0 and 
J = 2 —>■ 1 maps by Kang et al., submitted) are spatially 
well-defined in the pixel array and separation of out¬ 
flows and non-outflow pixels is easily performed. For this 
case the pixels that do not have any influence from the 
protostellar object, and which are included in the non- 
outflows combination, are Al, Bl, Cl, C5, D4, and E3 
of the SLW pixel array. Note that in some cases where a 
CO J = 1 —> 0 contour map was unavailable, the SPIRE 
field of view exceeded that of our comparison CO map; 
for these cases, there is some possible contamination by 
larger scale outflows. The line flu xes were t hen compared 
to the models produced bv lPon et all (2012:) to determine 
if the CO lines do follow the decay of turbulence trend, 
and to which properties each is most closely correlated. 

3.3. Grand Average Spectra 

3 http://irsa.ipac.caltech.edu/applications/DUST/ 

4 http://www.astro.princeton.edu/ draine/dust/dustmix.html 


When we compared the observational data to the mod¬ 
els, it was clear that none of the objects follow the PDR 
model. The data are a good match to the IPon et ahl 
(2012:) shock models, specifically to those of the highest 
(10 3 5 cm' 3 ) density, highest velocity (v = 3 km s _1 ) 
and higher magnetic field strength (B = 24 pG). How¬ 
ever, the only field that shows a clear difference between 
the outflow and the non-outflow pixels is around L1157. 
This could be because the object is oriented such that 
the pixels that contain outflows are clearly distinguish¬ 
able from those that do not (see Figure [2]); we may still 
be suffering from confusion with outflows in other fields. 
There could be additional influence by the outflows in 
the pixels which have been chosen as non-outflow pix¬ 
els and this may account for the close correlation of the 
fluxes in the outflow and non-outflow pixels around all 
but LI 157. In order to improve the signal-to-noise ra¬ 
tio, we computed a grand average of all seven objects to 
search for general trends. Spectra for all center pixels, 
outflows, and non-outflows were (separate ly) averaged 
together and compared with the I Pon et akl ( 2012 ) mod¬ 
els (Figure 0]). The CO J = 8 —>■ 7 through J = 4 —> 3 
line fluxes derived from pixels classified as outflow, non- 
outflow, and central protostar are indicated by purple, 
green, and light blue points respectively. The grand aver¬ 
age spectra for the outflow and non-outflow positions are 
shown in Figure [5] and Figure [6] respectively. The grand 
average spectra show clear distinctions between outflow 
and non-outflow pixels once the confusing line from [C I] 
at 370 pm has been removed, as seen in Figure 1. 



Wavelength (/jm) 

Fig. 5.— The flat spectrum averaged over all pixels encompassing 
the protostellar outflows for all seven sources in this study. 

These grand average line intensities are fitted by the 
highest density, hig hest v elo city and higher magnetic 
field model from IPon et all ( 201 2). The center pixels, 
outflows, and non-outflows all follow a distinctly differ¬ 
ent trend and the grand average of the non-outflow pix¬ 
els follows the shock models very closely. (See Figure E]). 
While there is still a concern about contamination from 
the outflows themselves, the different patterns for the 
distribution of emission over rotational levels suggests 
that the non-outflow positions are tracing dissipation of 
low velocity shocks rather than the high-velocity shocks 
in the outflow and center positions. If so, the fact that 
a model with higher density, higher velocity, and higher 
magnetic field compared to the isolated molecular cloud 
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Fig. 6. — The flat spectrum averaged over all pixels not contam¬ 
inated by protostellar outflows for all seven sources in this study. 


results of lPon et al.1 (12014 l2015bl iail in gas of similar den¬ 
sity, our data provide strong evidence for the dissipation 
of turbulence in molecular clouds. While substantial un¬ 
certainties remain in the timescale for dissipation, the 
data support simulations in which MHD turbulence de¬ 
cays within a crossing time. Particularly for a cloud like 
Taurus, with relatively distributed and low star forma¬ 
tion activity, the question of what inputs balance the 
decay presents an interesting challenge to theorists. 

The measured intensities for the grand average spectra 
of embedded protostars show clear differences between 
the average of positions in the outflow and positions not 
in the outflow. These results suggest that we are not 
simply seeing contamination from the outflow shocks. In¬ 
stead, we tentatively identify the CO spectra toward the 
non-outflow positions as a good match to models of tur¬ 
bulent decay in dense gas. Since these spectra are all in 
dense cores surrounding forming stars, the fact that the 
best fitting models have higher densities is not surprising. 


position fits these data suggests that the infalling en¬ 
velopes can be stirred by infall or by propagation of the 
outflow energy via magnetic fields. However, because 
the grand average spectrum is a mixture of data from 
different clouds, calculation of a quantitative dissipation 
timescale would be inappropriate. 

4. CONCLUSIONS 

We observe evidence of turbulent shocks in a region 
of the Taurus cloud far from protostell ar so u rces. This 
result supports the conclusions of iPon et al.1 (120141 . us¬ 
ing a different technique in a different cloud. Together 
with similar conclu sions in regions of more di ffuse gas 
(jCIodard et al.l[20T4 l and (jFalgarone et, al.ll2005H and the 
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